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\0 ' Abstract 

O' 



We determine the scaling equation of state of the three-dimensional 0{N) univer- 
sality class, for A^ = 5, 6, 32, 64. The N = 5 model is relevant for the S0{5) 
C^ ■ theory of high- Tc superconductivity, while the A^ = 6 model is relevant for the chi- 

, ral phase transition in two-color QCD with two flavors. We first obtain the critical 

Q I exponents and the small-field, high-temperature, expansion of the effective poten- 

r^ • tial (Helmholtz free energy) by analyzing the available perturbative series, in both 

fixed-dimension and e-expansion schemes. Then, we determine the critical equation 
k>( \ of state by using a systematic approximation scheme, based on polynomial rep- 

^ • resentations valid in the whole critical region, which satisfy the known analytical 



properties of the equation of state, take into account the Goldstone singularities at 
the coexistence curve and match the small-field, high-temperature, expansion of the 
effective potential. This allows us also to determine several universal amplitude ra- 
tios. We also compare our approximate solutions with those obtained in the large- A^ 
expansion, up to order 1/A^, finding good agreement for N > 32. 

Key words: 0{N) models. Equation of state. Universal amplitude ratios. 

Two-color QCD, harge-N expansion 

PACS: 75.10.Hk, 05.70.Jk, 05.70.Ce, 64.60.Cn 



Email addresses: agostino.butti@niib.infn.it (Agostino Butti), 
f.parisentoldin@sns.it (Francesco Parisen Toldin). 



Preprint submitted to Elsevier Science 



1 Introduction 



According to tlie Renormalization Group tlieory, it is possible to classify con- 
tinuous phase transitions into universality classes, determined by few global 
properties, such as space dimensionality, number of components of the order 
parameter, the range of interaction, symmetry group and symmetry-breaking 
pattern [1,2]. Some physical quantities, called universal, are independent from 
the microscopic structure of the interaction: models belonging to the same 
universality class have the same critical exponents, scaling functions and uni- 
versal amplitude ratios. One of the most important universality classes are 
the 0{N) classes, which are characterized by a A^-component order parame- 
ter, and a symmetry group 0{N) which is spontaneously broken to a subgroup 
0{N—1) in the low-temperature phase. For a recent review on this subject, see 
Ref . [3] . This universality class can describe the critical behavior of many phys- 
ical systems, which undergo a second order phase transition. For instance, the 
liquid-vapor transition in simple fluids and the magnetic transition in uniaxial 
(anti-)ferromagnets belong to the Ising universality class {N = 1). The XY 
model, which corresponds to A^ = 2, describes the helium superfluid transition 
and the Meissner transition in type-II superconductors. If the order parameter 
has A^ = 3 components, we get the Heisenberg universality class, which de- 
scribes the critical behavior of isotropic magnets. The three-dimensional 0(4) 
model is related to the finite-temperature phase transition in QCD with two 
hght fiavors [4,5,6,7]. 

In this article, we study the critical equation of state and universal quantities 
for the 0{N) three-dimensional universality class, with A^ > 4. We first con- 
sider the N = 5 and the A^ = 6 models. The 3-D 0(5) model is relevant for 
the so-called 5*0(5) theory of high- Tc superconductivity [8], since, according to 
this theory, the 5*0(5) symmetry should be approximately realized at a mul- 
ticritical point of the phase diagram ^ . According to universality arguments, 
the A^ = 6 case should describe the chiral phase transition in QCD with two 
fiavors and two colors [10]. This subject has been investigated by using an 
effective field theory in Ref. [11]. For all these 0{N) models, we first analyze 
the perturbative series available, in both fixed- dimension scheme [12,13,14] 
and e-expansion [15,16,17,18,19]; by resumming these series we determine the 
critical exponents and the coefficients rg and rg of the small-magnetization 
expansion of the critical equation of state in the high-temperature phase. In 
fixed-dimension we resum the series for the critical exponents ri{g), rjt{g) and 
for the beta function l3{g) up to six loops, and the series for r^lg) and rs{g) up 



^ However, it has been shown [9] that the 0(5) multicritical fixed point is in fact 
unstable: this implies that the 5*0(5) symmetry is not realized asymptotically at the 
multicritical point, but it controls the critical behavior in a preasymptotic region 
for t « 10-2. 



to four and three loops respectively. In e-expansion we resum the series for rj^e) 
and ri2{e) up to five loops and the series for r^i^e), rs{e) and for the universal 
quantity g4,{e) up to three loops. Then, we implement a systematic approxima- 
tion scheme for the equation of state, which provides an analytic continuation 
to the low-temperature phase. This procedure, based on polynomial represen- 
tations, has already been used for other models in the 0{N) universality class 
[20,21,22,23,24]. We determine, using this approximation scheme, the scaling 
functions and other universal quantities. We also compare the results with 
those in the large- iV limit, up to order 1/A^, that we calculate starting from 
the large- iV expansion of the equation of state reported in Ref . [25] . A recent 
review on the subject of large- iV field theory can be found in Ref. [26]. In order 
to check the convergence of this expansion, we repeat our work for two other 
0{N) models, N = 32 and A^ = 64. In these cases we find a good agreement 
of our results with those obtained in the 1/N expansion. 

This paper is organized as follows. In section 2 we summarize the scaling 
properties of the equation of state for the 0{N) model. 

In section 3 we introduce the notation and define some universal amplitude 
ratios that we calculate. 

In section 4 we introduce a parametric representation of the equation of state, 
valid in the whole critical region. 

In section 5 we consider a systematic approximation scheme, based on poly- 
nomials, that we use in our calculations. 

In section 6 we show critical exponents, universal couplings and other vari- 
ous universal quantities for the 0(5) and 0(6) three-dimensional models. We 
present a comparison of this results with those available in literature and 
with the 1/N expansion. We also plot the universal scaling functions we have 
obtained. 

In section 7 we consider two models with large A^: the A^ = 32 and the A^ = 64 
three-dimensional models. We report various universal quantities for these 
models, comparing the results with those obtained from the 1/A^ expansion of 
the equation of state. 

In appendix A we report the algorithm and formulas that we have used for 
the analysis of the perturbative series. 

In appendix B we report the 1/A^ expansion of the equation of state and the 
relevant formulas for the universal quantities. 



2 The scaling equation of state 



The equation of state is a relation between the magnetization M, the reduced 
temperature t = {T — T(.)/T(. and the external magnetic field H. It is obtained 
by noting that H is the derivative of the effective potential (Helmholtz free 
energy) A(t, M) with respect to M . Near the critical point H = 0, T = Tc, it 
can be written [1,2] in the scaling form 

H={B')-^MM^-^f{x), x = B^'HM-^"^, (1) 

— * 

where M = \M\, f{x) is a universal scaling function, fixed by the normal- 
ization /(O) = 1, /(—I) = 0, B'^ and B are the non-universal magnetization 
amplitudes at the critical isotherm and at the coexistence curve: 

M = B^HH^^^~\ t = 0, (2) 

M = B{-tf, t < 0, H = 0. (3) 

In the high-temperature phase t > 0, the critical equation of state can also be 
written as: 

H = ab^t^^F{z), z = hMr^, (4) 

where F{z) is a universal scaling function. Since the free energy is analytic in 
the plane (t, H) outside the coexistence curve and the critical point, F{z) can 
be expanded around z = in odd powers. The non-universal constants a and 
6, which appear in (4), fix the normalization on F{z) such that 

In the following, we refer to the universal coefficients r2„ as the universal 
couplings. Inverting (1), the equation of state can also be written as 

M = B,HH^/^-^E{y), y = {B / B,f'HH-^''^'^^\ (6) 

where E{y) is a universal scaling function. The relation between the scaling 
functions E{y) and f{x) is: 

E{y) = f{xr^/\ i/ = x/(x)-i/(^+^). (7) 

From (5) it follows the large-x behavior of f{x): 



f(x) =x''Yl fnX 



oo — 2n/3 
/ 
n>0 



The fact that the equation of state is analytical for iJ > implies that f{x) 
can be expanded around a; = 

oo 

fix) = 1 + E /n^"' (9) 

n=l 

and that for large- 2;: 

F(z) = / E Fr^-"/^. (10) 

n>0 

The relation between F{z) and f{x) is: 

z-'F{z) = F^f\x), z = zox-f, (11) 

where 

zo = hB = {Rtf' (12) 

is a universal amplitude ratio (see section 3). 

At the coexistence curve, that is for x -^ —1, spontaneous breaking of the 
0{N) symmetry occurs, so /(—I) = 0. In three dimensions, the presence of 
Goldstone modes leads to the prediction [27,25,28]: 

f{x)c^Cf{x+lf, x^-1. (13) 

This behavior is exact in the large-N limit [25]. The correction to (13) is less 
clear. If we define v = x + 1 and w = HM~^, using e-expansion it has been 
conjectured [28] that v admits a double expansion in powers oiw and ■u;(^~2)/2. 

V = cM'^-^^^^ + biw + C2w"'^^ + b2W^ + ... (14) 

In c? = 3 eq. (14) implies that v has an expansion in powers of w^^'^. However, 
this is not true in the large- A^ limit. In fact the 1/A^ correction to the leading 
behavior [29] shows the presence of logarithms for x ^ —1. 

From (13) one can derive that the transverse and longitudinal susceptibilities, 
near the critical point, for t < and iJ — > 0, diverge as: 

M , , 

XT = -^ (15) 

On 

In particular, in three dimensions xl diverges as H~^/^ for H —* 0. The scaling 
behavior of the longitudinal susceptibility in the critical domain can be derived 
from (6): 

XL = B,H'/'-'Diy), (17) 

where the scaling function D{y) is: 



Diy) ^ ] 



E{y) - JE'iy) 



Pf{x[ 



,1-1/5 



(36f{x)-xf'{xy 



This function has a maximum aX y = ymax > 0, which defines the crossover or 
pseudo-critical fine [3] (see section 3). 



3 Universal amplitude ratios 

Aside from scaling functions defined in the previous paragraph, there are other 
important universal quantities: the amplitude ratios. They can be defined 
from the critical behavior of zero-momentum quantities, such as specific heat, 
susceptibility, etc.. . . 

Near the critical point, the specific heat has the following singular behavior: 

T,CH = A^\t\-^ + b, (19) 

where the amplitudes A+ and A^ are for T > T^. and T < T^., respectively, 
and 6 is a no n- universal constant (we use the notations of Ref. [3]). The sus- 
ceptibility, in the high-temperature phase, diverges for t -^ 0+ as: 

X = NC+t-"'. (20) 

The critical behavior of the four-point connected function, at zero momentum, 
is given, in the high-temperature phase, by: 



m+^cn-'-'f'' . 



X4 = ^ 3 ' Ctt'^~'^\ (21) 

The divergence of the correlation length at the critical point is described by: 

e = f^t-^- (22) 

Universal amplitudes can also be obtained by considering the crossover, or 
pseudo-critical, line [3], which is defined as the reduced temperature tmax{H) 
where the longitudinal susceptibility has a maximum at H fixed: 

tmax{J^ ) ^ -'p-W , l2oj 

XLytmax: J^ ) = ^p'^max- 1^4) 

Corresponding to this maximum, the scaling variables x, y and z have the 
universal values Xmax, Umax, and Zmax, respectively. From equations (2), (3) and 
(19-23), one can define several universal amplitude ratios. Those involving the 
correlation length amplitude /''" can be obtained by using the zero-momentum 
four-point coupling (74, defined as: 



3iV X4 
N + 2x^e 



9^ - -TT—.Z^r (25) 



Table 1 

Definition of universal amplitude ratios and universal quantity (74 [3] . 



Uq = A+/A- 


Ra = {l- Uo)/a 


R+ = aA+C+/B^ 


Rt ^ -CtB^/{C+f 


R^ = C+B^-^/B^ 


R+ ^ (aA+y/^f+ 


P„ ^ T^B/B- 


P, ^ -T^^'C+/Ct 


Rp = C+/Cp 


9A ^ -Ct/[{C+?{f+f] 



At the critical point, g^ approaches a universal value, that we indicate with the 
same symbol. Several definitions of universal amplitude ratios [3] are reported 
in table 1, as well as the universal quantity (74. 



4 Parametric representation of the equation of state 



In order to study the equation of state in the whole critical regime, it is quite 

useful to introduce a parametric representation that implements the known 

analytical and scaling properties [30]. We introduce two variables, R and 9, 

such that: 

M = moRf^m{9) 

t = R{l-e^) (26) 

H = hoR^^h{e), 

where ho and mo are two normalization constants, i? is a nonnegative variable 
which measures the distance from the critical point. The functions h{6) and 
m{6) are odd and are conventionally normalized so that h{6) = 6 + 0{6^) 
and m{9) = 6 + 0{6^), for ^ — > 0. From (4), (5) and this normalization, one 
can obtain that ho /mo = ab'^; in order to take account of this, we define a 
parameter p such that 

ho = pab, mo = p/b. (27) 

As can be seen from (26), the line ^ = corresponds to the high-temperature 
phase, while on the critical isotherm 6=1. The coexistence curve is given by 
9 = 9o, the first positive zero of h{9). Clearly, consistence of the representation 
requires that m{9) > ior < 9 < 9o and h{9) > ioi < 9 < 9o. The 
mapping (26) is not invertible when its Jacobian vanishes; this happens when: 

Y{9) = {1- 9^)m'{9) + 2f39m{9) = 0. (28) 

Defining 9i as the smallest zero of Y{9), it must be 9i > 9o. 

All the scaling functions, scaling variables and universal amplitude ratios def- 
initions, that we have reported in sections 2 and 3, can be written in terms of 



the representation (26); the formulas can be found in [23]. We only note that, 
in order to reproduce the expected behavior (13) of f{x) near a; = — 1, the 
zero 9q of h{9) must be double: 

hi9) r^iO- 9oY, 9 ^ ^0- (29) 



5 Approximate equation of state 



Following [21], we now introduce an approximate polynomial representation, 
based on the parametric representation (26), that has all the expected proper- 
ties we enumerated in section 4. It extends the one considered in [20] and was 
already used in other three-dimensional 0{N) systems: the XY model [21,22] 
(iV = 2), the Heisenberg model [23] (A^ = 3) and the 0(4) model [24]. In this 
approximation scheme^ , we take h{9) as a 5th order polynomial and m{9) as 
a polynomial of degree 2n + 1: 

m{9) = 9(l+j2c^9A, 



m = 0[l-^. 



Z\^ (30) 

It can be shown that, in the N ^ oo limit, the scheme n = is exact [25]. 
The parameters p, 9q and n coefficients q, i.e. n + 2 parameters, are fixed 
by imposing that the equation of state, in the limit z — > 0, or, equivalently, 
^ ^ 0, reproduces the expansion (5). It is easy to see that we can implement 
the scheme at the order n if we know n+ 1 coefficients r2k'- this occurs because 
imposing the small-z expansion (5) at the ffist order gives a trivial equation. 

We see that there are two sources of error in this calculus. The ffist one is the 
uncertainty on the input parameters. The second one is the systematic error 
due to this approximation scheme. Indeed, if we are able to compute at least 
the n = and n = 1 schemes, we can estimate this error by looking at the 
differences between the results for n = and those for n = 1. Taking into 
account equations (30), we see that consistence of the whole calculus requires 
the coefficients Cj to be small. In particular, for the n = 1 case, it must be: 

]2 



cil^o < 1- (31) 



^ In [21,22,23,24] it was also considered another approximation scheme, denoted by 
(B). However, in this work, the scheme {B) didn't work, either by providing 9i < 6q 
or by failing to satisfy eq. (31). We mention that also in the A^ = 4 model, using the 
approximation schemes {B) and that of eq. (30), the latter provided better results 

[24]. 



6 Critical equation of state for the 0(5) and 0(6) models 



We present in table 2 the results for critical exponents and universal cou- 
plings for the models A^ = 5 and A^ = 6, obtained by an analysis of the 
fixed-dimension perturbative series [12,13,14] and of the e-expansion series 
[15,16,17,18,19]. In fixed- dimension we resum the series for the critical expo- 
nents ri{g), r]t{g) and for the beta function P{g) up to six loops, and the series 
for TqIq) and rs{g) up to four and three loops respectively. In e-expansion we 
resum the series for ri{e) and i]2{e) up to five loops and the series for re(e), 
rs{e) and (74 (e) up to three loops. We also show a comparison with other deter- 
minations available in literature and with the large- iV expansion. The series 
for rg, in the fixed-dimension scheme, is known up to three loops [14]: for this 
reason we can only provide a rough estimate of this coupling. The resumma- 
tion of universal couplings in e-expansion is constrained with the exact results 
available in dimensions d = and d = 1. The resummation is performed with 
the conformal mapping method; see appendix A for details. 

As explained in section 5, in order to compute the scaling equation of state 
and universal amplitude ratios, we need the input values of the critical expo- 
nents and universal couplings. We chose those final quantities by comparing 
the values obtained by resumming the fixed-dimension and the e-expansion 
series. For the 0(5) model our input values are r] = 0.031(3), z/ = 0.764(4), 
c/4 = 15.74(2), re = 1.71(2), rg = -0.3(5). For the 0(6) model they are 
7] = 0.029(3), u = 0.789(5), c/4 = 14.43(4), rg = 1.65(3), rg = -0.7(4). Since 
we have rg and rg, we can implement the approximation scheme up to the 
step n = 1. The results for universal amplitude ratios and other universal 
quantities are presented in table 3, along with a comparison with the 1/A^ ex- 
pansion. The errors reported are due to the variation of the input parameters 
in their error intervals. Inspecting table 3, we see that all the results obtained 
in the n = 1 scheme are consistent, within the errors, with the results for the 
n = case. In particular, for the A^ = 5 model, the central values for the 
n = 1 scheme are very close to those for n = 0: this is due to the fact that 
the predicted rg in the n = case matches the input value used in the n = 1 
scheme. Moreover, the condition (31) is satisfied. The equation of state for the 
0(6) model has been obtained by means of Monte Carlo simulation in Ref. 
[34]. The comparison should be done with some cautions, since the critical 
exponents in Ref. [34] are different from the ones we used in this work (see 
table 2). We report the results: R^ = 1.09(1), Umax = 1.34(5). Moreover, from 
the fits of Ref. [34] other universal quantities can be obtained: Rf = 7.6(2), 
F^ = 0.0233(9), Cf = 2.1(2), P„, = 1.13(2), P^ = 0.36(6). Comparing with 
our results shown in table 3, we see that there is a slight difference among the 
universal amplitude ratios involving the high-temperature phase. The scaling 
functions f{x), F{z), E{y) and D{y) are displayed, respectively, in figures 1, 
2, 3 and 4. In figures 3 and 4 we also show a comparison with the Monte Carlo 



Table 2 

Critical exponents, universal couplings and the universal quantity (74 for the models 
N = 5 and A^ = 6 obtained by resummation of the fixed-dimension series [12,13,14] 
(first row) and e-expansion series [15,16,17,18,19] (second row). A comparison with 
other determinations is shown. 



N 


54 


V 


z^ 


re 


^'s 


5 fixed-dim 


15.74(2) 


0.030(1) 


0.764(2) 


1.72(2) 


-1(3) 


e-exp 


15.6(1) 


0.034(2) 


0.764(6) 


1.70(1) 


-0.3(5) 




"15.69 


"0.034 


"0.766 


H.73 








^0.054 


^=0.784 


^2.88 


^^-13 


6 fixed-dim 


14.43(4) 


0.028(1) 


0.788(2) 


1.66(2) 


-1(3) 


e-exp 


14.23(8) 


0.032(2) 


0.790(6) 


1.64(2) 


-0.7(4) 




''14.39 


"0.031 


"0.790 


H.67 


^-0.28 




^14.59lfi 


"0.032(10) 


"0.804(3) 








'^umtu 


^0.034(9) 
»0.033(9) 
^0.037(9) 
^0.031(5) 
^.039(3) 


■^0.821(3) 
30.796(3) 
^0.819(3) 
*0.786(5) 
■'0.818(5) 










^^0.045 


^=0.820 


^=2.54 


fc_ll 



"[13] FT d = 3 exp, ''[14] FT d = 3 exp 

"[31] HT sc, ^^[31] HT bcc 

"[32] HTE sc unbiased, ■''[32] HTE sc 6'-biased 

s[32] HTE bcc unbiased, ^[32] HTE bcc 6'-biased 

*[33] MC, ^'[34] MC 

'^Large-A^ expansion to order 0{1/N), see appendix B 



fits of Ref. [34]. As can be seen, the curves n = and n = 1 are very similar. 
We observe a little difference between our curves and the Monte Carlo fits, in 
particular for large negative y in E{y) and for large positive y in D{y). This, 
again, can be interpreted as coming from the different critical exponents used. 
In fact, for y — > —00, E{y) ~ {—y)^: the discrepancies at large negative y are 
caused by the P = 0.425(2) exponent used in Ref. [34], obtained by their MC 
simulations, which is larger than our (3 = 0.406(3), obtained by field-theory 
methods. For y -^ +00, D{y) ~ R^y~'^. While our universal amplitude ratio 
R^ is compatible with that of Ref. [34], the critical exponent 7 = 1.555(10) 
used in this work is smaller than the one used in Ref. [34], 7 = 1.604(6). Fi- 
nally, we mention that the coefficient F^ appears to be very precise. This is 
not trivial, since the series (5) has a finite radius of convergence: starting from 
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Table 3 

Universal amplitude ratios and other quantities for the 



models N = 5 and A^ = 6. 



A^ = 5 


A^ = 6 




n=0 


n=l 


1/N 


n=0 


n=l 


1/N 


P 


2.16(1) 


2.16(9) 


- 


2.15(2) 


2.12(7) 


- 


01 


2.75(6) 


2.7(2) 


- 


2.56(7) 


2.5(1) 


- 


Cl 


- 


0.0006(234) 


- 


- 


0.008(19) 


- 


'•s 


-0.3(1) 


-0.3(5)* 


-13.46 


-0.5(1) 


-0.7(4)* 


-11.22 


no 


2.7(8) 


3(8) 


281.4 


5(1) 


8(8) 


234.5 


Uo 


2.21(3) 


2.2(2) 


1.47 


2.39(5) 


2.5(2) 


1.47 


Ra 


4.1(1) 


4.2(6) 


0.77 


3.8(2) 


4.0(5) 


0.72 


Rt 


0.279(9) 


0.28(2) 


0.595 


0.33(1) 


0.33(2) 


0.595 


Rt 


8.25(9) 


8.3(5) 


1.18 


8.3(1) 


8.4(4) 


2.99 


^x 


1.20(2) 


1.2(1) 


0.62 


1.12(3) 


1.15(9) 


0.681 


R- 


0.527(5) 


0.527(6) 


0.6185 


0.575(6) 


0.575(6) 


0.6573 


7700 

^0 


0.0214(5) 


0.0214(7) 


0.01961 


0.0196(5) 


0.0195(6) 


0.01750 


/J 


1.48(2) 


1.48(9) 


1.782 


1.57(2) 


1.55(6) 


1.819 


fo 


0.33(2) 


0.33(6) 


0.620 


0.41(2) 


0.40(5) 


0.684 


f^ 


-0.064(2) 


-0.06(2) 


-0.129 


-0.073(2) 


-0.07(1) 


-0.107 


c/ 


3.9(4) 


4(3) 


1.3 


2.7(2) 


3(1) 


1.3 


Pn. 


1.153(6) 


1.15(3) 


1.029 


1.133(6) 


1.14(2) 


1.040 


Pc 


0.330(2) 


0.330(4) 


0.4088 


0.319(3) 


0.318(5) 


0.3741 


Rp 


2.069(4) 


2.069(9) 


1.6593 


2.088(5) 


2.091(9) 


1.7827 


^max 


1.311(3) 


1.311(6) 


1.0194 


1.326(4) 


1.327(7) 


1.1077 


Xmax 


7.3(2) 


7.3(6) 


5.22 


6.7(2) 


6.9(4) 


5.18 


Umax 


1.44(2) 


1.4(1) 


1.089 


1.36(2) 


1.38(7) 


1.107 



* Input values 



the small- 2; (high-temperature) expansion, the approximation scheme provides 
an analytic continuation to the low-temperature phase. On the other hand, 
we see that the 1/A^ expansion provides values that differ significantly from 
the results obtained. 
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N=5 



N=6 





Fig. 1. The scaling function f{x) for the models A^ = 5 and A^ = 6. 
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Fig. 2. The scaling function F{z) for the models N = 5 and A^ = 6. 
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Fig. 3. The scaling function E{y) for the models N = 5 and A^ = 6. For the 0(6) 
model we show, for comparison, the results of Ref. [34] . 

7 Critical exponents and universal couplings for large A^ 



We present in table 4 the results for critical exponents and universal cou- 
plings for the models A^ = 32 and A^ = 64, obtained by an analysis of the 
fixed-dimension perturbative series [12,13,14] and of the e-expansion series 
[15,16,17,18,19]. In fixed- dimension we resum the series for the critical expo- 
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N=5 



N=6 




Fig. 4. The scaling function D{y) for the models N = 5 and A^ = 6. For the 0(6) 
model we show, for comparison, the results of Ref. [34] . 

Table 4 

Critical exponents, universal couplings and the universal quantity 54 for the mod- 
els A^ = 32 and A^ = 64 obtained by resummation of the fixed-dimension series 
[12,13,14] (first row) and e-expansion series [15,16,17,18,19] (second row). A com- 
parison with other determinations is shown. 



N 


54 


V 


V 


r% 


^8 


32 fixed-dim 


4.2295(5) 


0.00901(1) 


0.9584(2) 


1.1106(7) 


-1.3(2) 


e-exp 


4.209(3) 


0.0088(2) 


0.972(3) 


1.112(6) 


-1.244(7) 




M.230 


''0.009 


"0.958 


H.lll 


^-1.23 




'=4.22(1) 






^=1. 110(9) 


" - 1.24(2) 






'^0.00844 


'^0.966 


'^1.152 


^^-2.1 


64 fixed-dim 


2.23177(4) 


0.00467(2) 


0.9803(2) 


0.9819(4) 


-0.86(9) 


e-exp 


2.2280(3) 


0.00404(5) 


0.982(6) 


0.9846(3) 


-0.785(2) 






^^0.00422 


"^0.9831 


"^0.9929 


'^ - 1.052 



'='[13] FT d = 3 exp, ''[14] FT d = 3 exp 
"=[19] FT e-exp constrained 
Large- A^ expansion to order 0{1/N), see appendix B 

nents rj{g), rjt{g) and for the beta function [3{g) up to six loops, and the series 
for r^{g) and r^{g) up to four and three loops respectively. In e-expansion we 
resum the series for 77(e) and 772 (e) up to five loops and the series for r6(e), r^{e) 
and 5'4(e) up to three loops. We also show a comparison with other determina- 
tions available in literature and with the large- A^ expansion. The resummation 
of universal couplings in e-expansion has been constrained with the exact re- 
sults available in dimensions c? = and d = 1. The resummation is performed 
with the conformal mapping method; see appendix A for details. 

When A^ is large enough, the coefficients of the series become small. Thus, we 
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expect that the errors provided by the resummation technique become smaller. 
However, as we state in appendix A, the errors obtained by the analysis of the 
series don't have a statistical meaning: they are the variation of the quantity at 
hand, when the resummation parameters are varied. We believe that the errors 
provided by our procedure are substantially underestimated when A^ is large. 
In order to get the input values for the approximate equation of state with 
reasonable errors, we must again compare the results obtained with the fixed- 
dimension with those obtained with the e-expansion. For the 0(32) model 
our input values are t] = 0.0089(3), u = 0.965(8), g^ = 4.22(1), rg = 1.112(6), 
rg = -1.24(2). For the 0(64) model they are i] = 0.0044(4), u = 0.981(4), 
c/4 = 2.230(2), re = 0.983(2), rg = -0.79(1). 

In table 5 we present the results for various universal quantities, for the models 
iV = 32 and A^ = 64, along with a comparison with the 1/A^ expansion. 
The errors reported are due to the variation of the input parameters in their 
error intervals. Inspecting table 5, we see that, for both models, all the values 
obtained in the n = scheme are consistent, within the errors, with those 
obtained in the scheme n = 1. In particular, in many cases the central values, 
in the two approximation steps, are very close. Regarding the 1/A^ expansion, 
we note that, except for the coefficients r2„, even in the A^ = 32 model there 
is a good agreement with our universal quantities and those coming from 
1/A^ expansion of the equation of state. The agreement is even better for the 
A^ = 64 model: in this case almost all the quantities calculated in the n = 1 
approximation scheme are compatible, within the errors, with those obtained 
in the 1/A^ expansion. 
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A Resummation of perturbative series 



We have calculated the critical exponents and the coefficients tq and rg for 
the 0{N) model in three dimensions using the perturbative series from field 
theory available in the literature; these quantities are the input parameters 
for the determination of the critical equation of state and of the universal 
amplitude ratios, as explained in section 5. 
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Table 5 

Universal amplitude ratios and other quantities for the models A^ = 32 and N = 64. 



Af = 32 


Af = 64 




n=0 


n=l 


1/N 


n=0 


n=l 


1/N 


P 


1.9(1) 


1.90(7) 


- 


1.91(10) 


1.76(5) 


- 


01 


1.51(9) 


1.51(1) 


- 


1.48(8) 


1.406(7) 


- 


Cl 


- 


-0.0001(209) 


- 


- 


0.02(2) 


- 


'•s 


-1.2(3) 


-1.24(2)* 


-2.104 


-0.7(1) 


-0.79(1)* 


-1.052 


no 


19(7) 


19(3) 


44 


10(3) 


15(2) 


22 


Uo 


1.5(3) 


1.5(5) 


1.47 


1.4(2) 


1.5(4) 


1.47 


Ra 


0.5(3) 


0.5(6) 


0.51 


0.4(2) 


0.6(4) 


0.49 


Rt 


0.55(3) 


0.55(7) 


0.595 


0.55(3) 


0.58(5) 


0.595 


Rt 


10.4(1) 


10.42(4) 


10.310 


11.10(7) 


11.14(3) 


11.155 


^x 


0.939(2) 


0.94(1) 


0.940 


0.965(1) 


0.968(4) 


0.9701 


R- 


1.11(1) 


1.11(4) 


1.148 


1.40(2) 


1.42(4) 


1.447 


7700 

^0 


0.0092(2) 


0.0092(2) 


0.00892 


0.0081(1) 


0.00805(8) 


0.007934 


/J 


1.97(1) 


1.97(2) 


1.966 


1.984(5) 


1.983(7) 


1.9830 


fo 


0.94(2) 


0.94(3) 


0.941 


0.971(9) 


0.97(1) 


0.970 


f^ 


-0.020(6) 


-0.020(7) 


-0.0201 


-0.011(3) 


-0.011(3) 


-0.0100 


c/ 


1.04(2) 


1.04(3) 


1.051 


1.020(7) 


1.02(1) 


1.026 


Pm 


1.0800(5) 


1.080(3) 


1.0822 


1.0853(2) 


1.0860(9) 


1.08710 


Pc 


0.233(3) 


0.233(1) 


0.2334 


0.218(1) 


0.2177(6) 


0.21717 


Rp 


2.295(9) 


2.295(4) 


2.2843 


2.340(5) 


2.343(2) 


2.3421 


^max 


1.475(5) 


1.475(2) 


1.4664 


1.507(3) 


1.509(1) 


1.5078 


Xmax 


5.00(5) 


5.00(7) 


5.034 


5.01(2) 


5.01(3) 


5.017 


Umax 


1.171(2) 


1.171(8) 


1.1764 


1.181(1) 


1.182(3) 


1.1844 


* Input values 













Since perturbative expansions are divergent (asymptotic), a resummation pro- 
cedure is needed to obtain reasonable numerical estimates. Consider a generic 
quantity S{g) that has a perturbative expansion 



S{9) = E 'kg' 



(A.l) 



in the renormalized coupling constant g, corresponding to the coupling term 
of the 0^-Hamiltonian we are considering. We use a normalization for g such 
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that the beta function may be written as f3{g) = —g + g'^ + 0{g^) [13]. The 
large order behavior of the coefficients is generally given by 

Skr^k\{-afk^[l + 0{k-^)\. (A.2) 

Borel summability for g > ^ (proved in fixed dimensions d < 4 [35] for 0{N) 
symmetric models) requires a > 0. The constant a does not depend on the 
specific observable, unlike the constant c. From instanton calculations (see rfs. 
[1,36,37,38]) one can determine 

a = 1.32996798/(8 + N), (A.3) 

for the 0{N) model in three dimensions. We then introduce the Borel-Leroy 
transform B{t) of S{g) 

POO 

S{g)= / dtt'e-'B{gt), (A.4) 

Jo 

where 6 is a resummation parameter. Its expansion is 

This expansion is now convergent in the disk of the complex plane: |t| < 1/a, 
where a is the constant (A.3) that characterizes the large order behavior of 
S{g). In fact the singularity oi B{t) that is nearest to the origin is in ts = —1/a. 
However, in order to reconstruct the function S{g) from the integral in (A.4), 
we need to perform an analytic continuation of the expansion (A. 5). For exam- 
ple one can use Fade approximant to the series (A. 5): this is the Fade-Borel 
method. We have used instead a more refined procedure, the conformal map- 
ping method, based on the knowledge of the constant a: following Ref. [39], 
we perform an Euler transformation 

y{t) = ^^"\ (A.6) 

and rewrite the Borel transform as 

B{t) = —^--^j:bky{tr, (A.7) 

where a is another resummation parameter and the coefficients bk are deter- 
mined by the requirement that the expansion of B{t) in eq. (A.7) in powers 
of t is equal to the series (A. 5). If all singularities of B{t) belong to the real 
interval ]— oo, — tc], the expansion converges everywhere in the complex t-plane 
except on the negative axis for t < —tc- A good numerical estimate of S{g) is 
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therefore given by 

yigty 



Sig) c^Y^bk r dtt'e-'-y^^/-^, (A.8 

t^o -^0 [1 - y{gt)] 



where now the sum goes until k = p, since obviously we know only a finite 
number of terms of the perturbative series. 

We have also analyzed the e-expansion of the critical exponents and of the 
coefficients tq and rg, using the same conformal mapping method for the re- 
summation of the series in e (where e = 4 — d) and setting then e = 1. The 
constant a, which characterizes the singularity of the Borel transform, is given 
by 

a = -^—, (A.9) 

8 + N' ^ ' 

for the e-expansions in the 0{N) model [36,37]. 

Before resumming any series S{g) (or iS'(e)), we divided it by an opportune 
power of g (or e), so that the new series starts with a nonzero constant. 

In the framework of fixed dimension expansion {d = 3), we numerically re- 
summed the series for: 

• P{g) up to six loops [12,13]; we determined the IR fixed point g*: P{g*) = 0. 
We recall that the universal quantity g^, defined in table 1 and reported in 
tables 2-4, is related to g* by: (74 = 48'n'g*/{N + 8), for the three-dimensional 
0{N) model. 

• ri{g) and rit{g) up to six loops [12,13], using the value oi g*. The other critical 
exponents were determined through the hyperscaling relations among them. 

• rQ{g) and rs{g) up to four and three loops respectively [14]. 

As to the e-expansion we resummed the series for: 

• ri{e) and ri2{e) up to five loops [15,16]. The other critical exponents were 
determined through the hyperscaling relations among them. 

• g^i^e), 'rg(e) and r^^e) up to three loops [17,18,19]. 

The values of 5'4(e) and r2j{e) are known exactly in dimensions d = and d = 1; 
therefore we have also performed for these quantities constrained analyses 
using the same method described in [17,18,19] and references therein, which 
should improve the final results under the assumption of a sufficiently large 
analytic domain in e for the quantity to be resummed. 

Suppose that the exact values i?ex(ei), • • • Rex{ek) of the generic quantity R{e) 



17 



are known for a set of dimensions ei, . . . e^. Then define 



Q{e) = E 



i=l 



Rp 



e,: 



e- ei 



n 



and 



and finally 



S{e) 



R{e) 



n 



k 
i=l 



e,; - e,- 



QW 



-'limp I, f J 



[Q(e) + s{e)] n 



e — e,- 



(A.IO) 



(A.ll) 



(A.12) 



The resummation procedure 
tained by computing -Rimp(e = 
the values -Rex(ei), • • • -Rex(efc) 
that the e series which gives 
original one. In fact one can 
cients of S{e) decrease with k. 
should be smaller. 



is applied to S{e) and the final estimate is ob- 
= 1). When the polynomial interpolation through 
is a good approximation for -R(e), one expects 
the deviation has smaller coefficients than the 
check that, in the cases considered, the coeffi- 
Consequently also the errors in the resummation 



Any resummed quantity depends on the arbitrary parameters a and h: Sp{a, b), 
where p is the number of known terms in the original series. The parameters 
a and b may be used to optimize the resummation, see e.g. refs [39,40,41]. We 
shall now describe the algorithm we used; the idea is to choose the values of a 
and b providing the fastest convergence of the results in the number of loops. 



We tried to minimize numerically the quantity 

/ Sp{a,b) 



Sp^i{a,b) 



1 



(A.13) 



in a suitable range of parameters: —3/2 < « < 2 (for A^ = 5, 6), — 2 < a < 5 
(for A^ = 32, 64) and < 6 < 20 (for all the A^ considered). Sp-i is the quantity 
resummed using only p— 1 terms of the original series. When we found only one 
minimum of (A.13), we considered the positions of the minimum (oopt, ^opt) as 
the best resummation parameters. However in many cases, we found a whole 
line of zeroes for Sp—Sp-i in the plane {a, b), that is a line of minima for (A.13). 
In those cases we chose the optimal parameters (ftopt, ^opt) by minimizing the 
quantity 

fSp-i{a,b) 



\Sp-2{a,b) 
along the line of zeroes for Sp — Sp^i. 



I 



(A.14) 



Once we have determined the best resummation parameters, we considered 
different sets of approximants for the quantity to be resummed [41], obtained 
by varying the parameters in an opportune range a G [ctopt — Aa , aopt + Aa] 



and b G [6opt — A6, 6opt + A6] . One can take as the final estimate the average of 
5*^(0;, b) with a and b in this range (we considered half-integers a and integers 
b, after approximating ctopt and 6opt to the nearest half-integer or integer, 
respectively), while the error is indicated by the variance of Sp{a,b). It is 
not so clear how to determine the width of the intervals Aa and Ab [41]: in 
general the average of Sp is stable, within the quoted errors, while the variance 
strongly depends on the choice one is making. In order to obtain reasonable 
errors, we compared our results with those obtained by other authors in the 
case A^ = 4, and determined Aa and Ab so as to reproduce their errors. We 
have chosen Aa = 2 and Ab = 3; moreover we quote the error in the final 
results as two standard deviations. In this way we obtain for g* with A^ = 4 
the value g* = 1.375(4), which agrees with the result quoted by R. Guida and 
J. Zinn- Justin in Ref. [40]: g* = 1.377(5). 

For series in fixed dimension depending on g* we summed the squared er- 
rors due to the uncertainties of a, b, and of g*; the latter was estimated by 
{Sp{aopt, bopt, g* + Ag*) - Sp{aopt, &opt, 9* - ^9*))/'2. 

This algorithm is ad hoc and in some way arbitrary, but provides estimates 
consistent with each other. Moreover our results are in substantial agreement 
with those of other authors in the case A^ = 4 [40]. Therefore we consider our 
final values and error bars reliable, although we should be cautious in giving 
them the standard statistical meaning. 



B Large A^ behavior of the equation of state 

We consider now the large- A^ behavior of the critical exponents and of the 
equation of state for the three dimensional 0{N) model; our aim is to com- 
pare the results from the 1/A^ expansion with the numerical estimates for 
critical quantities obtained through resummation procedures and through the 
parametric approximation of the equation of state. 

The longest 1/A^ series for critical exponents available in literature are [42,43]: 
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16 1408 1024 



Nti^ iV297r4 Nm-K^ 



27 2, . 81 ,„ /1\ 61 2 683 



6 

(B.4) 



512 512 

+ 



A^Stt^ N^27ti^ iV3277r6 



9 ,2, o , 2r,./l\ 61 2 797 



2 ^ 8 ^ V2y 24 1^ 



32 32 /112 2 



(B.5) 
TT^ ) , (B.6) 



iV37r2 iV%4 V27 

where %l){x) is the logarithmic derivative of the F function. All the equalities 
(B.1-B.6) are understood to be true up to higher orders in 1/N . Since, in order 
to compute the large- iV expansion of the various universal quantities, we used 
the equation of state up to the order 1/iV, we have, for consistence, limited 
ourself to the order 1/A^ in the critical exponents. The large-iV limit of the 
equation of state and the 1/A^ corrections were calculated in Ref. [25]; since we 
will use them, we report their results in the case of three spatial dimensions 



f{x) = {x + lf + ^{x + l) 



' / X / X / X 12 (a;+ 1) log(a; + I] 
g{x) -{x + l) g{{)) ^ '-^ '- 



(B.7) 
where g[x) is given by 

1 r- 1 -\-2'K{\^x)^{k^)^^^ 



(B.8) 
and the functions i{k,r) and j{k,r) are defined as in Ref. [25] 



V, N 1 f k \ 
t[k,rj = — - arctan — 



j{k,r) 



Ank 

-3v^ 



(B.9) 



;7r (A;2 + r) (A;2 + 4r)' 



The function f{x) in (B.7) has the right normalization: /(—I) = and 
/(O) = 1. Starting from these equations, we derived the large-A^ expansion 
of universal amplitude ratios and of other quantities. 

The functions f{x) and g{x) are regular in a; = 0; therefore it is easy to 
compute the coefficients ff 

/N 2 + ^ (-^ - 2 ^(0) + 2 ^'(0)) , (B.IO) 

/2° = 1 + ^ (-^ - 2^(0) + 2^'(0) + /(O)) , (B.ll) 
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N [ 7r2 



/^-^(~ + «"(0) + ^ . (B.12) 



The values of g{x) and its derivatives in x = can be obtained from the 
definition (B.8) by calculating numerically the integrals in k. In the opposite 
limit a; -^ oo, we wrote the behavior of g{x) as 

8 8 Tc r°° 

g{x) = -{x + l) + —\og{x+l) + —\og-+ dkR{k,x), (B.13) 

TT^ TT"^ 4 JO 

where the function R{k, x) is 



R{k,x) = ———— + 



2(l + x) -l-27r{l + x)t{k,l) + 'jgffl 



(l + k)n^ (1 + P)^ ^2 (_L_ + (1 + p) ,(^^ 1)) 

(B.14) 
The last term in (B.13) can be expanded in powers of l/{x + 1) 



rdkR{k,x)=go + ^ + ^^^ + -^--^+o(-^—]. (B. 
Jo x + 1 [x + 1)^ (x + 1)'^ \[x + lpj 
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The coefficients gt can be obtained by Taylor expanding the integrand and 
calculating numerically the integrals in k. Then one can find the coefficients 

/r 

/o°^ = 1-^(1+^(0)), (B.16) 

/r^2 + lf-2 + .o-^-2.(0) + ii^l, (B.17) 



N \ "" TT 



2 -V ^ ' ^2 



/r = ^- (B.19) 

Using these results and the definitions of universal amplitude ratios, one can 
find: 



12 /-12+^o^+ 8 log(f ) 



i^l = 12 + ^( '^^"^°\""^^"^^M , (B.20) 



^x = 1 + 1^ (1 + ^(0)), (B.21) 

poo 11 TT^ (1 - ^0 + ^(0)) + 4 (3 + log(192) - 2 log(7r)) 

^° =144 + iV 72^^^ ' (^-^^^ 
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5 1 

6^N' 


5(2 + 5i7r2) 
37r2 






1 


35 5-2 








N 


6 


) 






1 


-35 


{2 + 6g27i^- 


-Q93 


2\ 
^ ) 



r, = 4^. (B.24) 



All these equalities are understood to be true up to order 1/A^. 

To derive the expressions of universal ratios involving the specific heat ampli- 
tudes, we need the effective potential (Helmholtz free energy), proportional to 
the universal function A{x). Using the notations of Ref. [3], we have 

f{x) = A{x) - ^A'{x). (B.26) 

(Ally 

A{x) is the solution of this differential equation which is analytic in x = 0. 
Solving (B.26) using the expressions (B.1-B.8) for f{x) and critical exponents 
up to 1/N terms, and then taking the limit A^ ^ 00, we obtained for A{x) 

7^2/24 \ 

A{x) = 1 + 3a; + 3a;2 - — + 3/(0) + g"\0) x\ (B.27) 

32 \ vr^ / 

Note that the knowledge of /(x) and z/ up to the order 1/N allows to determine 
only the leading term of the 1/A^ expansion of A{x). This is related to the fact 
that, in dimension d = 3, the critical exponent a tends to an integer for 
N ^ 00. Therefore, in the expression for the free energy in the large A^ limit, 
one cannot separate the term which is singular near the critical point from the 
analytic background. Ref. [26] gives the large- A^ behavior of the free energy 
and of the amplitudes R^ and Rt for generic fixed dimension d; as explained, 
these results cannot be extended to rf = 3. 

Using the expression (B.27), it is easy to obtain the large- A^ limit of the 
following universal amplitude ratios 



-;^(-^ + 3/(0)+/'(0) 
f/o = 32^V n ._ 3 28) 

^ ^ ' +3/(0)+/'(0) 



2 



32 V TT 

TT^ / 24 

Rt = -^ [--, + 3/(0) + g"'{0) ] . (B.29) 
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The leading behavior of Rt is 



R^ 

4:7!-' 



Rt={-^N] , (B.30) 



where R^ is given by its large- A^ limit (B.29). The 1/N series for the coefficient 
Cf is determined by the behavior of f{x) and g{x) near x = —1. We found the 
following expression 

g{x) = —{x + 1) log(a; + 1) + C{x + 1) + o{x + 1), (B.31) 

for a; ^ — 1. The constant C may be expressed as 
4 (2 + log(64) - 3 log(7r)) 



C 



7r2 



(B.32) 



^^ io [l + k " (1 + F) (4-l-P) arctan(|)j 
Then the coefficient C/ is (up to order 1/A^) 

Cf = l + ^iC-giO)). (B.33) 

The position of the crossover hne can be determined by imposing D\y) = 0; 
solving this equation up to order 1/A^, we found 



T 



=5 + — (64 - 30 ^'(5) vr^ + 300 g"{5) vr^ + 5 vr^ ^(5)) , (B.34) 

Vma. =— + -^ (320 - 150 g'{5) n^ + 1500 g"{5) tt^ + 60 ir^ g{0) 

65 A'157r"'65 ^ 

(B.35) 

+ 15 7r^^(5) + 16 log(216)), 



'-^yhj^ n ("^^^ + 5 ^0 ^^ + 30 ^'(5) n' - 300 ^(5) n' 



-5 7r^^(5) + 40 log(57r/4)j. (B.36) 

The universal amplitude ratios involving the crossover line are 

^™ =7? + ..on L^ . (320 - 150 g'{5) n' + 1500 g"{5) n' + 60 vr^ g{0) 
65 A/ 30 v5 65 iT^ ^ 

+ 15 TT^ ^(5) + 528 log(6) - 600 log(5)) , (B.37) 
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p, 



3125 3125 



-100 + 12 TT^ + 3 c/o TT^ + 30 g'{h) tt^ - 300 /(5) tt^ 



" 15552 A^466567r2 

-3 7r2^(5) - 192 log(2) - 144 log(3) + 168 log(5) + 24 log(7r^ 



(B.38) 



Rv = Y + j^ (-76 + 15 TT^ + 15 g\h) vr^ - 300 /(5) vr^ + ISO log(5/6)) 



We report the results after calculating numerically the integrals: 



^6 



rio 



6 ' N 


-67.314 


N 


1406.83 



N 



Uo = 1.4674 



R^ = 0.4674 



R+ = 0.594715 



RJ = 12 



^x = l 



54.0813 

N 

1.91184 9.22 



A^ A^2 



i?^ = 0.361715 ATs 



7-1 OO 

^0 



1 0.0633049 



144 



N 



/? = 2 



/2° = 1 



/3° 



c/ 



1.08789 

N 

1.89798 

N 

0.642744 

N 

1.64657 



P™ = 1.092 - 



0.313824 



Pr = 0.200939 



N 
1.0391 



A^ 



Rr, = 2.4 - 



3.70369 



N 



1.54919 



2.64874 



N 



Xri 



5 + 



1.08647 



A^ 



1.19247 



0.515247 

TV ■ 



In the last equations we quoted from [44] the correction l/N"^ to the universal 
amplitude ratio R^. 
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